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Abstract 

QQ ' We discuss a quantum effect in tlie diffusion process by developing a theory, which takes the finite 

■ curvature of the potential field into account. The transport coefficients of our theory satisfy the 

Tlj" , well-known fluctuation-dissipation theorem in the limit of Markovian approximation in the cases 



of diffusion in a fiat potential and in a potential well. For the diffusion along a potential barrier, 
the diffusion coefficient can be related to the friction coefficient by an analytic continuation of the 
fluctuation-dissipation theorem for the case of diffusion along a potential well in the asymptotic 
time, but contains strong non-Markovian effects at short times. By applying our theory to the 
^ ' case of realistic values of the temperature, the barrier curvature, and the friction coefficient, we 



show that the quantum effects will play significant roles in describing the synthesis of superheavy 
elements, i.e., the evolution from the fusion barrier to the conditional saddle, in terms of a diffusion 



X 

d ' process. We especially point out the importance of the memory effect, which increases at lower 

temperatures. It makes the net quantum effects enhance the probability of crossing the conditional 
saddle. 
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I. INTRODUCTION 

Diffusion processes take place in a variety of problems jl[ Q, 14 ) 1^ ls| • The concept has 
recently been applied to theoretically describing the synthesis of superheavy elements 0, ■ 
It is now well accepted that it is not sufficient for the two nuclei in heavy-ion collisions 
to overcome the fusion (i.e., the Coulomb) barrier in order to form a heavy compound 
nucleus such as superheavy elements. Since the conditional saddle is located inside the 
fusion barrier for heavy-ion collisions between two heavy nuclei, two nuclei have to further 
progress inwards to approach inside the conditional saddle. The idea of Refs. 0,0, called 
fluctuation-dissipation dynamics, is to describe the time evolution from the fusion barrier 
to the conditional saddle as a diffusion process. 

Though this approach is very attractive and is offering much useful information, one 
needs to examine the applicability of one of the basic assumptions made so far, i.e., the use 
of the standard fluctuation-dissipation theorem which holds at high temperatures to relate 
the diffusion coefficients to the friction coefficients. Since superheavy elements are stabilized 
by shell correction energies, one has to synthesize them at reasonably low energies, as low 
as 1 MeV or below. On the other hand, the curvature of the conditional saddle is also of the 
order of 1 MeV. It is thus required to carefully study quantum effects. An interesting issue 
is to explore the connection between the diffusion and friction coefficients in the diffusion 
process along a potential barrier under such circumstances. Though the generalization of the 
Einstein relation to the case of diffusion along a potential well is well known, the modification 
in the case of a diffusion along a potential barrier has so far been discussed only in a limited 
number of papers 0, 0, 1^ . 

The aim of this paper is to examine this quantum effect caused by the finite curvature of 
the potential barrier and by the low temperature aspect of the diffusion process. To this end, 
we first develop a novel quantum diffusion theory which leads to a Fokker-Planck equation 
with non-Markovian transport coefficients. We then apply it to the situation relevant to 
the synthesis of superheavy elements. We will show that the quantum effects, especially the 
non-Markovian effects, are very important. 

In Sec. II, we briefly sketch the derivation of the Fokker-Planck equation with non- 
Markovian transport coefficients, which include quantum effects. In Sec. Ill we discuss the 
quantum effects on the diffusion coefficient. In Sec. IV we apply the formalism to analyze 
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the quantum effects by assuming the parameters which are relevant to describe the diffusion 
process in the synthesis of superheavy elements. We summarize the paper in Sec. V. 



II. QUANTUM DIFFUSION EQUATION WITH NON-MARKOVIAN TRANS- 
PORT COEFFICIENTS 

We derive the required diffusion equation by extending the quasilinear response theory 
developed in Ref. |lO| so as to include the curvature of the potential barrier. Details of the 
derivation will be published elsewhere Here we sketch the main steps: (1) We start form 
the von Neumann equation for the system consisting of space A for macroscopic degrees of 
freedom and space B, for microscopic degrees of freedom. (2) We introduce the classical 
trajectory given by {q{t), pit)), t being the time and q and p being the coordinate and the 
conjugate momentum operators of the macroscopic degrees of freedom. (3) We move to 
the Galilei transformed coordinate system specified by q{t) and p{t). (4) We keep only up 
to the second order terms of the potential and the coupling Hamiltonian in the expansion 
with respect to the fluctuations around the classical trajectory. (5) We solve the coupled 
equations describing the A and B spaces in the lowest order approximation concerning the 
fluctuating force. (6) We make a Wigner transform of the resultant extended von Neumann 
equation for subspace A. 

Denoting the Wigner transform of the density operator of the subspace A by Daw, we 
finally obtain 

—DAw{p,q,'t) = -Y7P»^ — ^C"?' 



i-E) d (_o) d (+o^ d 




dpadqp 



The non-Markovian property of the diffusion process due to quantum effects, more specifi- 
cally the effects of the curvature of the potential barrier, is hidden in the transport coeffi- 
cients, which are generalized from the 1st and 0th moments of the response and correlation 
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functions, "* ^"^^ X^^\ 6.g. as 



X%''\t) = f dt,S{t,t,)xi-,\t,t,), (2) 



'to 

xiy\t) = f dt,C{tM)x^:^{tM)- (3) 

J to 

The C and S are defined by using the curvature of the potential C as 

C(t,ti) = cos[fi(t-ti)], (4) 

s{t,u) = -^sin[^](^-^l)], (5) 

with Vt = ^JcJM when C > 0, i.e., for the motion in a potential well, and 

C(t,ti) = cosh[fi(t-ti)], (6) 
S[tM) = -^^=sinh[fi(t-ti)], (7) 



with = a/|C|/M when C < 0, i.e. for a motion along a potential barrier. The response 

-)(t,ti)and 



and the correlation functions, XaJi^^^^) ^^"^ X^aJit^ti)^ are given by 



xi~J{t,h) = TTrBi[Ut)Jf^ih)]DBih)), (8) 



X%\t,h) = ^Trs([/„(t),/^(ti)]+DB(ti)), (9) 
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in terms of the fluctuation force fait) defined by 

fait) =uUt,to)FaUB{t,to), (10) 

with 

In Eq. (fTT|) . pcit) is the density operator of the total system in the Galilei transformed 
coordinate system. The ^5(^,^0) is the time evolution operator of the subspace B, which 
satisfies the partial differential equation 

d 

ih—UBit,to) = hBit)uBit,to), (12) 
with the initial condition UBito,to) = 1. In Eq. (I12j) . the effective Hamiltonian is given by 

hBit) = HBix) + VMt),x), (13) 



where Hb{x) is the unperturbed Hamiltonian of space B and Vc the couphng Hamiltonian. 
The density operator D{t) describing our basic equations [(H)), (jH)) and Q] is defined from 
the density operator in the Gahlei transformed coordinate space pcit) by 

PG{t) =UB{t,to)D{t)u^Bit,to). (14) 

The Dsit) is defined by Dsit) = TrA{D{t)). 



III. QUANTUM EFFECTS ON THE FLUCTUATION-DISSIPATION THEOREM 



The time evolution of the subspace B should in principle be determined by solving Eq. ()12p 



or the corresponding reduced von Neumann equation lol 
density operator by the canonical distribution, 

PBit) ^ e^^{p{t)[J^ -liBit)]}. 



12j |. Here we approximate the 



(15) 



We further replace /^^(t) by Hb in Eq. ()15|1 to be consistent with the linear response the- 
ory. One can then easily show by introducing spectral function ^, that the following 
well-known generalized fiuctuation-dissipation theorem follows in the limit of Markovian 
approximation for the motion along a potential well, 



A+E) 



it) 



A-o) 

X_(y.a 



with the effective temperature given by 



it) 



MT\ 



(16) 



-Ml coth 



'-P{t)hVL 



T (T > hO) 
\Ml (T < Ml). 



(17) 
(18) 



Equation (jT)) is then nothing but the Kramers diffusion equation [l^ postulated in Refs. 0, 
Is I, though there exist some modifications such as the temperature being replaced by the 
effective temperature. 

As declared in the introduction, our interest in connection with the synthesis of super- 
heavy elements is the diffusion process along a potential barrier instead of along a potential 
well. In this case, one needs to specify a model in order to further discuss the properties 
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of the diffusion coefficients. We assume the Feynmann- Vernon model Q], which has been 
used also by Caldeira and Leggett ,1^] to discuss macroscopic quantum tunneling, assume 
the Ohmic dissipation, and use the Drude regularization [3| by introducing the following 
cutoff function for the spectral density of the environment, i.e., the sub-space B: 

, . 1 



1 + 



We define 



Y{t) 



X 



i+E) 



(t) 1 



(19) 



(-0) M' 

Xoo 



(20) 



where Xoo is the expression which the odd-moment of the response function takes in the 
limit of the Markovian approximation or in the asymptotic time. The Y{t) consists of three 
terms, two of which strongly depend on time. 



_:^^2g-c(t-to)cot 

4 



X 



1 



Pit) 



n=l,2,. 



X 



7r2?i 

hl3{t) 



7r2n 



00% 



(21) 



1 




(22) 



The third one, which gives the asymptotic value, reads 



0^2 



Ml cot 



'-j3{t)hVt 



(23) 



Note that the right-hand side (r.h.s.) of Eq. is the analytic continuation of the r.h.s. 
of Eq. (fT7|) concerning the sign of the barrier curvature if one ignores the minor change due 
to the first factor. Such analytic continuation formula of the r.h.s. of Eq. (1171) has been 
argued in Ref. p. Our theory contains additionally non-Markovian effects. We note that 
Y{t) reduces to the temperature at very high temperatures. 

Figure HI shows the time dependence of Y{t), i.e., the ratio of the diffusion to the asymp- 
totic friction coefficients, for four different values of the temperature. The barrier curvature 
VL and the cutoff frequency Uc in the spectral density in the Drude regularization have been 
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FIG. 1: Ratio of the diffusion to friction coefficients. 



fixed to be HQ = 1.0 MeV and hw^ = 20.0 MeV. The horizontal hnes show the position 
of each temperature. The figure shows that the classical fluctuation-dissipation theorem 
postulated in the original Kramers paper ^ and that also in Refs. 0, 0] holds only at very 
high temperatures. The more striking result is the non-Markovian effect, which appears as 
the strong variation of Y{t) with time. We remark that the span where Y{t) shows this 
strong variation gets longer roughly in proportional to the inverse of the temperature as 
lone as T < huJr- 



IV. BARRIER CROSSING PROBABILITY 

We now apply our formalism to discuss quantum effects in the diffusion process from the 
fusion barrier to the conditional saddle in the synthesis of superheavy elements. 

We first note that the average values of q and p are zero in the Galilei transformed space, 
and that the solution of Eq. is a Gaussian. Therefore, one can set 



DAw{p,q,'t) 



1 



27rA 



-exp 



(^qq'^pp '^qp^ 



(24) 
(25) 



where yi = q and 1/2 = P and the 2x2 matrix a is the inverse matrix of the 2x2 matrix 



^qq ^qp 



^qp ^pp 



(26) 



which determines the fluctuations, i.e., the mean square deviations from the average values. 
The values of aij are obtained by solving the following coupled equations given by Eq. (jT]): 



d 
dt 



I 



_2_ 

M 



-X 



i-O) 



J_ 

M 



\ 



+ 



( 








(-0) 



\ 



(27) 



The Wigner distribution function for the macroscopic motion in the original space fixed 
frame is given by 



PAw{q,P,t) = DAw{q - qit) , P - Pit) , t) , 



(28) 



once the Wigner distribution function in the Galilei transformed space Daw is obtained 
We represent the conditional saddle by a parabola, 

1 



Ks(g) 



(29) 



and calculate the probability to cross the conditional saddle in order to form a compound 
nucleus by 

1 ^ [q-qit)r 



Pit) 



dq- 

^y2'^■aqq(t) 



exp 



2 I v/2M^) 



(30) 



We ignore the radial dependence of the friction tensor. Denoting the initial position and 
momentum of the classical trajectory of the macroscopic variable as (go, Po), the position at 
time t is given by 



qit) 



e 2 



qo 



(cosh —t + sinh —t] +2^ sinh —t 



(31) 
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FIG. 2: Comparison of the probability to cross the conditional saddle calculated by quantum 
diffusion equation and by assuming the classical fluctuation-dissipation theorem. 



with j3' = \/ (3'^ + 4^2. We adopt the value of the reduced friction parameter (3 from previous 
studies of fluctuation-dissipation dynamics using Langevin equation as /5 = 5 x lO^^s^^. 
The curvature of the potential barrier is assumed to be HQ = 1 MeV, which is relevant to 
heavy nuclei. The initial position go is chosen to make the height of the conditional saddle be 
4.0 MeV, and the mass parameter to correspond to the reduced mass in the colhsion of the 
mass number 48 and 238 nuclei. We defer the study of the effects of purely non-Markovian 
terms ^ind x*^^*^-* and leave them in determining the <Jqq{t) in the following analyses. 

Figure 121 shows the probability to cross the conditional saddle as a function of the initial 
kinetic energy which is measured relative to the height of the conditional saddle V b- We 
remark that the ratio q{t)/ ^J2aqq{t) in Eq. (jHIH) converges to an asymptotic value. It was 
used to evaluate the probability to cross the conditional saddle P. We choose three values 
for the temperature. The solid lines are the results of our theory, while the dot-dashed 
lines are the results when the classical fluctuation-dissipation theorem has been assumed 
by ignoring the quantum effects due to the finite curvature of the conditional saddle. The 
figure clearly shows that the quantum effect is important at low temperatures, which are 
relevant to the synthesis of superheavy elements. 

Our theory contains a memory effect. In order to discuss the connection to a previous 
work jo^, we artificially isolate the memory effect by calculating the probability to cross the 
conditional saddle by using the asymptotic value of the diffusion coefficient. The result is 
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FIG. 3: Analysis of the memory effects. 



added in Fig. |21by the dotted hne. We observe that the probabiUty to cross the conditional 
saddle is reduced by the quantum effect if one ignores the memory effect. This is because the 
asymptotic diffusion coefficient in the quantum theory is smaller than that obtained from 
the classical fluctuation-dissipation theorem. A similar effect has been shown in Ref. jsl. 
Our study shows in addition that the memory effect overcomes this effect and finally the 
net quantum effects enhance the transmission probability of the conditional saddle. In other 
words, the net quantum effects reduce the fusion hindrance. 

In passing, we wish to mention that quantum effects on diffusion process are also discussed 
in Ref. following a different approach. However, there are some important differences 
in the expressions of transport coefficients. For example, the diffusion coefficient given by 
Eq.(8) in Ref. does not seem to match with our asymptotic formula (23), as well as those 
presented in Refs. ^, 0] . We also wish to refer to Ref. [l^ , which discusses the dynamics of 
barrier penetration in a thermal medium for the inverted harmonic oscillator by using the 
influence functional formalism of the path integral method. 



V. SUMMARY AND FUTURE DEVELOPMENTS 



We have presented a diffusion theory which takes the finite curvature of the potential 
field into account. The theory is then applied to the case where the temperature, barrier 
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curvature, and the friction coefficient are taken to represent a realistic situation of the diffu- 
sion process from the fusion barrier to the conditional saddle in the synthesis of superheavy 
elements. We have thus shown that the quantum effects will play an important role. We 
especially pointed out the importance of the memory effect which has been omitted in any 
previous works. It makes the net quantum effects enhance the probability of crossing the 
conditional saddle, while the quantum effect reduces it if the memory effect is ignored, as 
has been shown in Ref. ^j. 

We have artificially left out some of the genuine non-Markovian terms, i.e., the odd 
moment of the correlation function and the even moment of the response function. Also, we 
have assumed a sharp distribution at the initial time and left out the effects of spreading 
of the initial distribution. We will discuss these effects as well as the dependence of the 
quantum effects on the strength of the dissipative force in forthcoming papers. One of 
the interesting problems is to clarify whether our conclusion concerning the role of non- 
Markovian effect is special to our specific choice of the Caldeira-Leggett model, especially 
to the Ohmic dissipation, or holds in general. This is another issue which we will explore in 
the near future. 
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